MCTP-08-57 



Holographic Geometric Entropy at Finite Temperature 
from Black Holes in Global Anti de Sitter Spaces 



Ibrahima Bah^, Leopoldo A. Pando Zayas^ 
and Cesar A. Terrero-Escalante^ 



^ Michigan Center for Theoretical Physics 
Randall Laboratory of Physics, The University of Michigan 
Ann Arbor, MI 48109-1040 



^ Departamento de Fisica, Centro de Investigacion y Estudios Avanzados del IPN, 
Apdo. Postal 14-740, 07000 Mexico D.F., Mexico 



Abstract 

Using a holographic proposal for the geometric entropy we study its behavior in the 
geometry of Schwarzschild black holes in global AdSp for p = 3, 4, 5. Holographically, 
the entropy is determined by a minimal surface. On the gravity side, due to the presence 
of a horizon on the background, generically there are two solutions to the surfaces 
determining the entanglement entropy. In the case of AdS^ , the calculation reproduces 
precisely the geometric entropy of an interval of length / in a two dimensional conformal 
field theory with periodic boundary conditions. We demonstrate that in the cases of 
AdSi and AdS^ the sign of the difference of the geometric entropies changes, signaling a 
transition. Euclideanization implies that various embedding of the holographic surface 
are possible. We study some of them and find that the transitions are ubiquitous. In 
particular, our analysis renders a very intricate phase space, showing, for some ranges 
of the temperature, up to three branches. We observe a remarkable universality in the 
type of results we obtain from AdS^ and AdS^. 



1 Introduction and outlook 



A very natural construction in field theory is the geometric entropy. Given a system in a 
state |\E' >, one defines a quantity associated with that pure state and a geometrical region 
A by forming the pure state density matrix p = >< \E'|, tracing over the field variables 
outside the region to create an impure density matrix pA = Tr^p, where B is the complement 
of A. One then evaluates the von Neumann entropy Sa = — Tr^ip^ log pa associated with 
this reduced density matrix. 

One of the first calculations of this quantity was performed in the eighties [1] for the case of 
a scalar field propagating in a black-hole background. This paper noticed its proportionality 
to the area of the black hole horizon. A similar problem but in a strictly field theory sense 
was posed by Srednicki. Namely, take a ground-state density matrix of a free massless field 
and trace over the degrees of freedom residing outside an imaginary sphere. The fact that 
the associated entropy is proportional to the area was connected to the physics of black holes 
[2] and considered as a concrete computation supporting the idea that the entropy area law 
is a more general phenomenon that originally anticipated. 

A lot of activity during the nineties was focused on precise computations of the geometric 
entropy in various situations [31 HI [5] . The fact that the geometric entropy came out to be 
proportional to the area elicited a lot of speculations about its role in the black hole entropy. 
Perhaps most of that activity was crystalized in a claim by Callan and Wilczek attempting to 
clarify the relationship between the geometric entropy as defined above and the Bekenstein- 
Hawking entropy of a black hole. Namely that geometric entropy is the first quantum 
correction to a thermodynamic entropy which reduces to the Bekenstein-Hawking entropy in 
the black hole context. More precisely, the Bekenstein-Hawking entropy and the geometrical 
entropy are the classical and first quantum contributions to a unified object measuring the 
response of the Euclidean path integral to the introduction of a conical singularity in the 
underlying geometry. Some concrete computations followed to support this claim, including 
the very similar statement: The appropriately defined geometric entropy of a free field in fiat 
space is just the quantum correction to the Bekenstein-Hawking entropy of Rindler space 

mm 

The geometric entropy has turned out to be an important measurement of the entan- 
glement between the two subsystems, it has been extensively studied in this context which 
includes its use as an order parameter for exotic phases (TU]. Recently, there has been con- 
siderable interest in formulating measures of quantum entanglement^ and applying them to 
extended quantum systems with many degrees of freedom, such as quantum spin chains. 

The AdS / CFT correspondence states roughly that some field theories are mathematically 
equivalent to string theories containing gravity [11]. The question of entanglement was 
addressed in this context by many authors, including [121 [13] • A crucial breakthrough was 
taken by Ryu and Takayanagi [HI [T5] who proposed a holographic way of computing it and 

^Entanglement entropy and geometric entropy have the same mathematical definition. The main dis- 
tinction is the situations in which they have been applied. In this note we use entanglement and geometric 
entropy interchangeably. However, we keep in mind that quantum transitions take place at zero temperature 
and therefore, in our context, we favor the term geometric entropy. 
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showed that it precisely matches expectations from two dimensional field theories. Many 
more works in this area has followed [IS] -(SO]- More recently in [HT], a prescription has been 
suggested for a holographic computation of the geometric entropy. The prescription for the 
holographic computation of Sa consists on evaluating the area of a minimal surface, in the 
dual supergravity background, whose boundary is the subsystem A: 



The entropy is obtained by minimizing the above action over all surfaces that approach the 
boundary of the subsystem A. This is a paradigm shift for no longer the geometric entropy 
is treated as correction to the black hole entropy as intuitively done in various papers from 
eighties and early nineties; the statement now is that the geometric field theory entropy can 
be calculated as a classical quantity in a supergravity background. 

Of particular interest is the case when the supergravity background has a horizon itself. 
Various groups considered this situation [12], [I3]. Recently, we revisited the question of 
entanglement entropy in the presence of horizons in which cases there are two surfaces 
that minimized the above action and therefore a competition ensues among them that we 
interpret as phase transitions in the field theory (see figured]). One of the surfaces represents 
a continuous configuration while the second one is piece-wise smooth and extents from the 
boundary to the horizon. In appendix |X] we prove that this last surface is also minimal. 

For a typical class of embedding and for many types of solutions we showed in [32] that 
generically there is no transition. The main result of this paper is that after Euclideaniza- 
tion, if we consider a slightly different class of embeddings, motivated largely by [31], then 
transitions are rather typical for AdS^ and AdS^. 



Figure 1: Two competing configurations for the entanglement entropy in the presence of a 
black hole horizon. The green surface represents a continuous configuration while the red 
surface is piece wise smooth and starts by going straight down from the boundary to the 
horizon. The subsystem A is given by the blue section. Its characteristic length is /. 
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1.1 Universality of transitions 

Let us state what we consider to be the main result of our paper. The phase diagram of Ad S 4 
and AdS^ is universal, the differences being in the quantitative values. We briefly summarize 
it here. The metric of the Schwarzschild black hole in global AdSn+i is given as: 

ds' = f^{r)dT^ + 4^ + r'dQl_, (1.2) 



/n(r) 



where 



M 

Mr) = ^-^.+T- (1-3) 

,^2 \de^ + sui^md^'^ + siYi^i^)d(l?) for n = 4 

\d9^ + sm\9)d(lf' for n = 3 ^ ^ 

In the case of AdS^, when we fix the 0-cycle we can compute the geometric entropy b 
embedding the surface as r(r) or r{9). In the case of AdS^, we consider r(r) and r(?/;) 
Since r, 9 {ip) have finite period, there exist a maximum / which we denote by Imax and 
represent with a dashed orange line in the plots below. 

For the r(r) embedding, we observe that there exist two critical values of the radius of 
the black hole horizon or, more physically, two critical temperatures, for which the behavior 
of the difference of the areas, AA{1) = Acontinuous — Apiece-wise, drastically changes. We 
tabulate these different behaviors below and denote yo = r horizon/ L the value of the radius 
of the horizon in units of the AdS radius. 

1. yo > ^o"*^; has one branch. In this case, it is smooth and negative as seen in figure 



2(a) So we do not observe any phase transition. 



ycriti > > yo'^^'^i has three branches and can change sign. In this case, has 
two points where it fails to be different iable. The first and the third branches intersect 



creating a loop as seen in figure 2(b) This suggests that AA jumps from branch one 
to branch 3. Thus we observe a first order phase transition. Furthermore, this occurs 
at a negative value of AA, which suggests that this transition occurs in the phase 
dominated by the smooth surface configuration. 

Vo < 1/0"*^ 5 AA has two branches, one positive and the other negative as seen in figure 



2(c), The branch where AA is positive occurs for / > Imax] it is not physical. The 
physical branch however, does not go to zero as / approaches Imax- This suggest that 
A A discontinuously jumps to at / = Imax- 

This phase structure is observed in various contexts and seems to be very universal. For 
example, it is found in the thermodynamics of charged AdS black holes [33], in the study 



^ Later on this paper we discuss this kind of cycles and their relationship with the more typically used 
temporal cycle. 
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(a) 2/0 > (b) yg"" > 2/0 > yS"" (c) Vo < 2/§"" 

Figure 2: For the embedding r(r), we plot against / at different ranges of Hq. The dashed 
orange hne is the maximum allowed value for /. 

of Nambu-Goto strings in gravitational backgrounds dual to field theory with flavors [34J- 
More importantly, this is the typical phase diagram of a van Der Waals ga^. 

When we consider the r{6) {r{ijj)) embedding, we observe a generic transition between 
the smooth and piece-wise smooth surfaces, since AA always changes sign at a value of / less 
than Imax] see figure [3l 
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Figure 3: AA is plotted against / for the embedding where the radius depends on an angu- 
lar direction. We generally observe transition between the smooth and piece-wise smooth 
surface. The dashed orange line is the maximum allowed value for /. 



1.2 Organization of the paper 

The paper is organized as follows. In section [2l we discussed the geometric entropy in AdS^ 
and find perfect agreement with the field theory calculation. Essentially, we holographically 
compute the geometric entropy of a segment of length / in a two dimensional conformal field 
theory with periodic boundary condition on a circle of radius Im- Sections [3] and H] discuss 
the holographic geometric entropy for the Schwarzschild black hole in global AdS coordi- 
nates. There are various possibilities for the embedding of the minimal surface depending 
on the finite temperature boundary conditions and the form of the dependence of the radial 

•^We thank C. Nunez from highlighting this univcrsahty to us. We hope to elaborate on this apparent 
universality in a joint collaboration. 
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coordinate; we explore a few of them. In section [5] we present some comments about the 
imphcations of choosing different cycles for the surfaces. We conclude in section El 



2 Geometric entropy in AdSa 

One of the most interesting results of [13] is the holographic calculation of the entangle- 
ment entropy for two dimensional conformal field theories. This interesting result has been 
expanded to a fairly elaborate choice of the subsystem A, namely, a discontinuous set of 
intervals in [22]. The starting geometry is the BTZ black hole 



ds^ = f(r)dT^ + -^dr^ + r^d(j)^ (2.1) 
fir) 

where 

f{r) = l-m + — = l-m + y^, (2.2) 

with y = r/L. In the Euclidean version of the metric we have the option of considering 
(j) as the Euclidean temporal direction. The subsystem A is then defined as the segment 
< r < where (3 is the period of the r cycle and / < 1. We consider minimal surfaces 
whose embedding is defined as r(r) with boundary conditions r(r — > 0) = r(r — > pi) — > oo. 
The metric of the induced surface is then, 

ds' = fir){l + ^)dT' (2.3) 

with area 

Ac= / vW)\n + ^dT. (2.4) 



The surface is minimal if r(r) satisfies 

fir) 



fir) fin) 



- 1. (2.5) 



In the above expression r^, is the turning point of the smooth surface corresponding to 
r'(r) = 0. The area can be conveniently rewritten as an integral over y 

Ac = / -j^dr = 2L / dy. (2.6) 

JO Vfir*) Jy, Vfiy) - fiv*) 

For the manifolds we are considering, there is another solution to the boundary problem 
which consists of a piece-wise smooth surface. The area of such surfaces is given by 

dr 

^p = 2/ (2.7) 
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The difference in area is then 



dy 



dy 



(2.8) 



The separation at the boundary / depends on the turning point of the surface in the bulk y^ 
as follows: 

l='Jir dy (2 9) 

These expressions can be integrated analytically to yield 

AA = 2L\n(^\ /(yj = -arcsin f^V (2.10) 

\y*J TT Vy*/ 

from which we obtain 

A^(/) = 2Lln(sin(7r/)). (2.11) 

Note that AA < 0. Therefore, we observe the entropy must be given by the smooth surface 
in AdSs. 

We can rewrite the expression for the entropy in terms of intrinsically field theoretic 
quantities. Recall that the central charge in the two dimensional CFT can be related to 
Newton's constant G and the radius of AdS3 (which is denoted by L) as: c = 

Since the holographic entanglement entropy is given as = we obtain 

= ^ = ^ ln(sin(vr/)) = ^ ln(sin(vr/)). (2.12) 

Interestingly, equation (12.101) shows that the length / is bounded by 1/2 and at that 
value the difference in entropies vanishes. The difference of entropies is symmetric under 
/ ^ l — l, this can be interpreted as meaning that the entropy of the subsystem A defined by 
the segment I and its complement, defined by the segment 1 — / are the same. The equality 
Sa = Sb, for the entropy of a subsystem A and its complement B, is a general property of 
the entanglement entropy and it is true in higher dimension. It is also worth noting that the 
difference in entropy is independent of yo, making it independent of the periodicity along 
the original temporal direction r. This property changes when we go to higher dimensions. 
In fact, it will allow for an interesting phase structure. 

Our holographic result perfectly agrees with the entanglement entropy of a subsystem of 
length / in a finite system with periodic boundary condition [TJ [32] and finite length. The 
precise field theoretic expression is 

5(0 = ^lnf^sin(f^))+ci (2.13) 

3 \7ia Irn J 

where Im is the size of the system, ci is a non-universal constant and a is the lattice spacing 
corresponding to a UV cut-off. 
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Our result from the continuous surface is 

A. - 2L ,n ^ S(l) ^ | In .„(.o) ^ In (1^ .n(.o) + f (A) . 

(2.14) 

This expression captures the field theory answer fl2.13p precisely. 



3 Geometric entropy from AdS4 

In this section we analyze the geometric entropy in the holographic context of the Schwarzschild 
black hole in AdS^^. The dual field theory is expected to be a 2 + 1 superconformal field 
theory at some inverse temperature /3 living on a the conformal boundary of AdS^ which is 
basically a 2-sphere. Although the precise conformal field theory can be understood in terms 
of the AdS/CFT correspondence, our main motivation comes from the fact that there seems 
to be a sort of universality for field theories in 2 + 1 dimensions that make some of our results 
relevant for a wide class of problems. For example, systems whose gravity dual asymptote 
to AdSi^ have been recently considered in connection with various condensed matter systems 
[3911101111]. 

The Euclidean geometry is given by: 

1 M 

ds^ = f{r)dT' + —-dr^ + r\d9^ + sm\9)d(P^), f{r) = 1 + — . (3.1) 

j[r) r L'^ 

Using the convenient coordinate y = r/L, we have: f{y) = ^{y^ + y — m) where m = M/L. 
Defining the position of the horizon as the largest root of f{yo) = we find that to avoid 
conical singularities the dimensionless periodicity in the r direction must be 



L l + Syi' 



/3=T = T^- (3-2) 



In the Euclideanized form of the metric we can clearly choose to embedded the surface 
at fixed r or at fixed ijj. We have considered embedding with fixed r previously [32]. In this 
paper we focused on embeddings where another coordinate is kept fixed. Namely, considered 
embeddings with ijj fixed. Then we have two options for the embedding of the surface, r{6) 
and r(r). We proceed to explore both of them in what follows (see Section O for comments 
on these embeddings). 

3.1 r{9) minimal surface 

The induced metric for this embedding is, 

ds' = f{r)dT^ + r\l + -^)de\ (3.3) 
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where the prime means derivative with respect to 6 and < 6 < ttI. The area is given by: 

tI 



Ac = (3 







The equation of motion can be written as 

r'^fijr) 



- 1. 



-de. 



(3.4) 



From the equation of motion, we can rewrite the area as, 



(3.5) 



/5 



n^/fin) 



d9 = 2(3 



dr. 



(3.6) 



The area of the piece-wise smooth surface that also solves the boundary condition problem 



IS 



Ap = 2p dr = 2[3{r^-ro). 



(3.7) 



The difference in areas, Ac — Ap^ is 



yvm 



Vy'fiy)-y*fiy*) 

The relation between the size of subsystem A, I, and the minimum of the smooth surface y^: 



IS 



Ky* 



y*Vfiy*) 



--dy. 



(3.9) 



'y* yV fiy)Vy^fiy) - ylfiv*) 

Our goal is to show that the difference in entropy changes sign. We managed to do that 
analytically by studying some limits of the above expression. First, we consider the limit 
when the smooth surface drops into the bulk almost all the way to the horizon. In this case 
we can show (see appendix [B]) that 



AA{y, = + e) > 



(3.10) 



where e « y^. In the other limit, when the surface remains in the large y region, we can 
simplify the expression using x = y/y^, to the form: 



2PLy, 



X 



x^ 



l\dx-l 



2(3Ly^ [0.4009 - 1] < 0. 



(3.11) 
(3.12) 



We thus observe that the difference in area changes sign. However, this is not enough to 
claim transition as the value of / must be bounded in order to have a well defined problem 
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on a circle of finite size. Given our normalization we can claim the existence of a transition 
only when there exist a yl such that AA{yl) = and the length / evaluated at that value 
of satisfies l{y^) < 1. We can then supplement the above statement by showing that 
the function /(y*) < 1 for all black hole sizes, that is, for all values of yo- This is observed 

0,0.1,0.5,1. In fact, it is numerically 
0, which corresponds to global 



in figure 4(a) where / is plotted against y^: at yo 
observed that 1 > /(y*, y^) > l{y^, yl) if y^ < y^. When yo 



AdS, there is no conical singularity problem; thus the period of r can take on any value. If 
the period is not zero, it is clear that AA will change sign since the above argument only 
requires a non-zero period for r. 





(a) 



(b) 



Figure 4: 4(a) shows the function l{y^:) at yo = 0,0.1,0.5,1. For the r{6) embedding, we 



observe that / is strictly bounded by 1. 4(b) shows the behavior of A A as a function of /. 
We generically observe a transition. 



We conclude that the difference in area changes sign for all finite yo- For completeness 



we also plot, in figure 4(b), the difference in areas as a function of the separation length /, 



that is, AA{1) for various values of yo. 

Let us finish this subsection with a comment about the ultraviolet behavior of the entropy. 
In the large limit, (13. 6p and (13. 9p integrate to 



A, = 2PL^ y, 



v^r(i 



-y* 



2 1 0Fr(: 
^y* r(i) 



The expression for the area of the smooth surface in this limit is 



(3.13) 



A,{1) = 2PL^ y. 



r(i) 



(3.14) 



We can then write the entropy, S 



AG 



(4y, in terms of the Planck constant, /p, and the 



number of M2 branes, iV, from the relations 



G 



(11) 

N 



(4) 



G 



(11) 

N 



RlM 



2L = R 



■s7 



7"7 



(3.15) 
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The entropy is then 



where C = ' 



. We notice that there is a hnear divergent piece, we might identify Voo/lp 

with an ultraviolet cutoff. It is interesting to note that the power of N'^^^ is different from 
the naively expected N^^"^. The simplest reason comes from dimensional analysis, that is, 
having the new scales (3 and Too in the problem makes us deviate from the conformal result of 
N'^^"^. Namely, each length scale should have contributed N^^^ in the conformal limit, we can 
simply verify that iV3/2 ^ ^7/6^^2/6^ g^^g optimistic that iV^/^ versus N^/'^ might have a 
real physical explanations in terms of degrees of freedom involved in the entanglement. In 
fact we will see similar deviations from the standard conformal scaling in all the problems 
that we tackle in this paper. 

The entropy from the piece-wise smooth surface is also 

where we have used the fact that ro = ^1 + ^1 — 4^7^ j ^ • 

We observe that the divergent piece is the same in both expected. Furthermore 

this quantity is proportional to the area of the boundary of subsystem A, (Sr^o- This is in 
agreement with the area law discussed in [15] and [2]. We see that Sp does not depend on /. 
Therefore, when the entropy difference changes sign experiences a discontinuity. 

3.2 r(r) minimal surface 

The induced metric for this surface is given by 

ds^ = /(r)(l + '^)dT^ + r^de\ (3.18) 

where the prime stands for derivative with respect to r and < r < The area of the 
embedded surface 

Ac = 71 I r^/fi^)Jl + j^dT. (3.19) 







The equation of motion for the minimal surface is then 

r'^(r) r'^fijr) 



1. (3.20) 



with area 



Ac = 71 ^j^^ dr = 271 dr. (3.21) 
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The area of the piece-wise smooth surface is also, 



Ap = 2tx 



rdr 



ro»L 



2'iTLrnr, — C 



(3.22) 



where C = 

Ap as 



1(1-/0 



00 



x-^/^-Vx3- i^^ ) ^ 0.144. We then obtain the difference in areas 



AA = 2txL^ 



y* 



The length / of the subsystem in this case is 



dy 



y* ydy 



yo 



(3.23) 



^ _ 2L rv- 



y*Vfiy*) 



-.dy. 



(3.24) 



P Jy, f{y)Vy'fiy)-y'J{y*) 

Let us first consider the limiting regions of the above expressions. We observe that 

AA{y, = yo) = 0. (3.25) 

For large y^,, we can write 



AA 



27rL^ 



y* 



X 



1 dx 



yo 



ydy 



2'KL'^y^ [0.4009 - 1] < 0. 



(3.26) 
(3.27) 



In the last line, we used the fact that the second integral is in the large limit. 

The analysis above does not show whether the difference in entropy changes sign nor does 
it show whether there exist a transition. In figure [5] and O we have plotted / vs. and AA 
vs. / at appropriate values of y^. We generically observe that /(yo) = \ from the numerical 
analysis. In appendix Owe proved analytically that as yo the length / 1/2. This 

feature is consistent with the fact that yo is really the origin of this background and since 
it is non thermal, the difference in area must satisfy AA{1) = AA{1 — I). Thus we should 
not consider cases for / > 0.5. This was also observed for AdSs. Furthermore, numerical 
analysis show interesting behavior for AA{1). We observe that AA is negative and monotonic 
for yo > 0.217 (see figure 6(c)), then from yo = 0.217 to yo = 0.189, AA and / both have 
three branches. Within this interval, AA can change sign. Then from yo = 0.189 to zero, 
AA and I both have two branches. In this interval, the branch where AA is positive exist 
only for / > 0.5. In figure U\ we show snapshots of AA{1) as yo decreases in the interval 
0.217 > yo > 0.189. We observe that AA has a first order phase transition without changing 
signs; then when yo < 0.189, it experiences a discontinuity from some negative value to 0. 
This phase structure is observed in the thermodynamics of charged AdS black holes |33j . 
fundamental strings in supergravity backgrounds dual to field theories with fundamental 
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(a) 



(b) 



Figure 5: In 5(a) we have plotted / against at various values of yo and in 5(b) , we have 



zoomed in the region around Z = .5. 





(a) 



(b) 



< -0 




(c) 



Figure 6: In 6(a), we plot the difference in area for values of yo less that .189. In 6(b), we 
have zoomed in the region around / = .5 to show interesting transition structure. In 6(b) 
we plot the difference in area for values of yo greater than .217. 
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Figure 7: We plot the difference in area for r(r) in AdS4, vs. / at difference values of yo. These 
snapshots track the the difference in entropy as it experiences a first order phase transition. 



flavor [M] and Van der Waals gas. In upcoming work in collaboration with C. Nunez we 
will elaborate on this analogy further. Here we conjecture that AA behaves as a Gibbs 
potential or free energy depending on which cycle we fix. In this picture, / behaves like 
inverse temperature. 

From equations (13.211) and (13 .24^ we find the following ultraviolet behavior for the area: 



^(0 



27rL' 



v^r( 



r(i) 



The entropy of the smooth surface for small / is then 



Sc{l) 



n7 (327r^iV)^/3 
(2^ Tp ' 



r(i) 



2(27r) 



The entropy from the piece-wise smooth surface at small (3 is 



S„ 



^7 (327r2A^)^/3 n-r {?,2Ti^Nf'% 



(27r 



(2vr) 



P 



(3.28) 
(3.29) 



(3.30) 



(3.31) 
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The ultraviolet behavior is linearly divergent, as in the previous embedding. The power 
of N'^/^ differs from the N"^/^ of the previous subsection and we note this as the result 
of breaking the conformal invariance. Namely, the presence of r^o explicitly indicates a 
missing factor of A^^/®, indeed, the conformal behavior would have been observed otherwise 
since N'^/"^ = N^^^N^^^. As expected, the divergent piece is proportional to the area of the 
codimension 3 boundary. This boundary has length 27rroo. To leading order, the finite piece 
is proportional to l/l. 



3.3 General Observations 

From the above analysis, we demonstrate that the holographic geometric entropy always 
indicates a transition for configurations in the Schwarzschild black hole in AdS^. We found 
that for the minimal surface parametrized by r{6) a transition occurs for i/q > 0, that is, 
for any nonzero temperature. We proved this by showing that there exist a y'^ where AA 
changes sign. It was then observed that at y'^ the value of the length is in the appropriate 
range: / < 1. This implies that there is a l"^, namely l'^ = l{yl)i where AA{1) changes sign. 

For the r(r) embedding, we observe a more interesting structure. The difference in area 
have three phases enumerated below. 

1. In the interval 0.217 < yo < oo, AA is negative, and / is bounded from above by 0.5. 
The result here is similar to the AdSs. 

2. In the interval 0.217 > yo > 0.189, AA and / has 3 branches. With this range, AA can 
change sign. 

3. In the interval 0.189 > yo > 0, AA and / has 2 branches and AA is positive only for 
/ > 0.5. This suggests that AA must jump from some negative value to zero at / = 0.5 
since the embedding exist only for / < 0.5. 

Although we have found the critical values of the dimensionless radius of the black hole 
horizon ?/q^**^ and y^"^"^ numerically, it is useful to determine them more systematically. The 
first critical value, y^"^^, occurs when I and AA are no longer monotonic. This implies these 
functions have local minima and maxima. Thus the critical value is the maximum of yo from 
which one starts to see these local extrema. At such value, the equation 

dl 

— = 3.32 

dy* 

must have a real solution. At yg^**^, l{y*) must also have an inflection point. We must also 
have 

So at the first critical value yo"^^, is given by the solution of the following equations 

dl d^l 

— = 0, and — = 0. 3.34 

dy* dyi 
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The second critical value, occurs when /(?/*) transitions from having 3 branches to 2. This 

dl 
dy* 

dl 



occurs when one of the zeros of is i/q"^'^. Thus at the second critical point, we have: 



= 0. (3.35) 



4 Geometric entropy from AdS^ 

Consider the Schwarzschild black hole in global AdS^ with the following metric 

ds^ = f(r)dT^ + -^dr^ + idO^ + sin\e)idi)^ + sin^r^)^^^)) (41) 
fij) 



with 



M 

/(O = 1 - H + Z^- (4.2) 



Introducing y = r / L and m = M/L'^, we have 

TD 

f{y) = 1- +y^. (4.3) 

We consider the entropy when (p is fixed. The periodicity in r is set to /5 to avoid conical 
singularities in that direction, its dimensionless periodicity is 

4.1 Minimal surface 

First we evaluate the smooth minimal surface. The induced metric is 

ds^ = f{r)dT^ + {r%'' + r^sm^{e))d^^ + r^d9^, with h= ""'^j^ , (4.5) 

f{r) 

The area is then 

^c = P J r^^W)\lh^ + sin2(^)c/^# = (5 j r^^/f^p{b)d^P, p{h) = j ^Jb^ + sm\e)de, 

(4.6) 

where p{b) is an elliptic function. The equation of motion is 

rl^^IKMO)=r'yWMb), with qib) = T ^^E==de. (4.7) 

The value of r^, corresponds to r' = or 6 = 0. After switching integration from ip to r, we 
obtain 

-p{b)dr. (4.8) 
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The piece-wise smooth surface evaluates to: 

A, = 2(3 rdrn = 7r/5(r^ - r^) (4.9) 

Jro 

The difference in areas is 

-{p{b) -TTb)dr + p7i{rl-rl). (4.10) 

Now we want to set up the problem so that we can integrate over b. The equation of motion 
is not so readily solvable to give r(6, r^) due to the presence of /(r). In the y coordinate, the 
equation of motion can be written as 

y' + y'-my'-ptg\b)=0 (4.11) 

where 

Pl = yl^/Iiy*), m = y',+yl and 9{b) = ^. (4.12) 
As integrals over b, I and AA are 

AA = 2L'pJ^^^y^y^ipib)-nb)db + L'Pniy',-y'^), (4.13) 

^ ^ 2 r dy{y,,b) db ^2 dyjb_ 

Trio by/y^ + y^ - m tt Jq ^ dbbg{b)' 

One solves for ?/(?/*, 6) from (14. lip and then evaluates these integrals. Furthermore, by using 
(14. lip , one can write 

dy ^ ptgjb) dg{b) 1 

db y db 3y' + 2y^ - yl - yi ^ ' ' 

from which one obtains, 



= 2L'P r P*f^ I , , (p(6) - .b)db + L^p.{yl - yl) 

Jo b db 3y^ + 2y^ - y^ - y^ 



2_pI ridgi^ db_ 

TT io b db 'iy^ + 2y^-yl-yl 



I = t!!± (4 16) 



As in the case case of AdS^, we explore two regimes determined by y*, that is by the turning 
point of the surface. For y^, yo, it can be shown that 

AA(|/o + e) >0. (4.17) 
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In the other hmit, that is, for large y^^ we have ~ and y^ ~ utd^^^i^)- This allows us 
to write AA, in the large y^, limit, as 

= 2L^^J^^y^^±{p{b)-^b)db-L^n(3y! (4.18) 



Since 



7o & db Sg^/'iby"^^ ' ' 2 



< (4.19) 



This shows that changes sign for some ranges of y^. Similar to r{9) in AdS^, this is not 
enough to conclude that the difference in area generically changes sign. We need to show 
that I is appropriately bounded. 

Let us present the ultraviolet behavior of the entropy. In the large limit, we can write, 

A = 2L'8^"^M^^p((,)*. y'«,.y"W (4.21) 
= PL'^Txyl,' L^0vil-iO7... (4,22) 
Using the expression for / in this limit 

2 dg 1 „ 0.416 



we obtain 



= (3^rl-L'(3^. (4.24) 



We can write the entanglement entropy, S = ^Tgy, in terms of field theory parameters from 



the relations 

gS^°) = 8A.V^ = Stt^^, L' = 4na"X, G'-^^ = (4.25) 

The entropy from the smooth surface for small / is 

_ {47r)'/'Q, 1 ^2 ^ _ (47r)VU 21.8(3 1 
' 327r5 /3/4 13 327r6 AV4 vr^/, Z^- ^^-^^^ 

where Ig is the string length scale Ig = ^/oH . We note a few interesting general properties that 
we reproduce. First, there is the quadratic divergence that is typical to four-dimensional field 
theories. For example, Srednicki had already discussed it in [2]. It was also more recently 
discussed and obtained for various field theories, including their holographic formulation by 
Ryu and Takayanagi in [15]. The 1//^ dependence of the finite piece is also rather general in 
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four dimensional field theories, as widely discussed in [15]. Both these statements seems to 
follow from dimensional analysis but it interesting that they are confirmed by the holographic 
calculations. As expected, we also obtain the A'^^ factor in the expression for the entropy. 
However, unexpectedly we have powers of the 't Hooft coupling /. We explained these powers 
as effects of breaking conformal invariance by introducing some scales in the problem, in this 
case (3. Recall that in the standard holographic calculation of the entanglement entropy 
there are no powers of A [T3] . We find this dependence on A very peculiar and it would 
be interesting to pursue its interpretation. We can also write the entropy coming from the 
piece-wise smooth surface. This entropy is independent of / and is given as 



327r^ 



13/4 ''c 



1287r4 (3 ' 



1 



2/?2 



1/2N 



(4.27) 



Here we have to replace rg for /3 through the relationship. 



2/5 



7r2L2 



(4.28) 



for a thermodynamically stable black hole. Notice that there is a maximum f3 for which the 
black hole exist. 





(a) 



(b) 



Figure 8: 8(a) shows the function /(?/=„) at yo = 10 ^, 0.1, 0.5. 8(b) shows the behavior of AA 



as a function of /. We generically observe a transition. 



We plot some results that verify numerically that we always get a transition. We note, 
as we did in the introduction and outlook section, that the result is qualitatively the same 
as for the r{6) embedding in AdS^. 
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4.1.1 r(r) Minimal surface 

The induced surface in the Euclidean frame is given as 

ds' = f{r){l + '^)dT^ + r^dnl. 
The subsystem is given by the segment < r < /?/, the area is given as 



A = 2ti 



Pir] 



dr. 



where prime is derivative with respect to r. The equation of motion is then 
The of the surface can then be written as 



A 



27r 



-.dr = 8 



TT 



'0 rl^/J{7~) 

The piece- wise smooth surface is given by r =constant, with area 

.2 / 



dr. 



Ap = in 



r 



-.dr 



'^'^Lr^ - TT 



In 



L2 + 2rl 



The difference in area is then 



Voc 



We can also write the relationship l{y^] 

, 2L 2 fjj-^ dy 

I = -p-y*vf{y*) 



dy 



yo 



dy 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



f{y)Vy'fiy)-ytfiy*y 

Let us now consider the limits of AA. Similar to the AdS4^ case, we immediately obtain 

AA{y, = yo) = 0. (4.36) 
In the case when y^, is large, we can rewrite the integral as 

,2 



AA 



AttL' 



yl 



X 



— X ] dx 



y 



yo 



dy 



0.28440y2 _ 



y, y2 



ya 



dy 



< 0. 



(4.37) 
(4.38) 
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In the first integral, we used the fact that for large f{y) ^ and then introduced the 
coordinate x = y/y^,. The inequality in the last line follows from the fact that for large y^ the 
second integral diverges as |?/^. The change in the difference of areas from zero to negative 
does not, by itself, indicates a transition. A similar situation was observed in the case of 
AdS^^. Although more work is still required at this point, we verify explicitly that the result 
for AdS^ is qualitatively the same as in 74^5*4, the difference is in the actual numerical values 
of the critical temperatures, that is, of the critical values of y^. In figure [HI we plot snapshots 
of ^A{1) in the range of y^ where we observe first order phase transitions. 

The ultraviolet behavior of the entropy can be obtained as before. In the large limit, 
we need the large y limit of /(?/), fiy) ~ y"^- We thus obtain 

A, = 27rL\^-27rL=^ln(2/oo) -47rLV0.216, / = J-0.431, (4.39) 

py* 

which translates into the following result 

A^{1) = 2nLrl-2nLHn[^) -AnL'^-^. (4.40) 

We can write the entanglement entropy, S = ^^sj? in the small I for the smooth surface as 



(4.41) 



SvrS /52 /2 

The entropy of the piece- wise smooth surface can be obtained from (14.331) . it is given by 



's 

, , w^awXy'VA , f, (, '213'' V"' 



In addition to reproducing the area law, we observe that the divergence here have a 
logarithmic term. 



5 Comments on the various cycles 

In the previous sections we have considered various configurations where we have fixed a 
periodic coordinate (j). In previous work [32] we considered fixing the temporal cycle which 
we denoted by r in this paper. A natural question then arises as to what is the relationship 
between these two cycles and what are the rules for choosing one or the other. 
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(a) 2/0 = .3 (b) zjo = .286 (c) yo = .266 




(g) yo = .212 (h) 2/0 = .19 (i) yo = .01 



Figure 9: We plot the difference in entropy for r(r) in AdS^ vs. / at difference values of 
Uq. These snapshots track the the difference in entropy as it experiences a first order phase 
transition. 

It is of course, a well known fact that the presence of two potential time cycles usu- 
ally induces transitions in field theory. Generically there is a low temperature phase and 
a high temperature phase. Perhaps the cleanest example is Witten's reinterpretation of 
the Hawking-Page phase transitions [121 There the competition between two saddle 
points of the action, namely, the Schwarzschild black hole in AdS and the thermal AdS 
was interpreted as the supergravity dual of the confinement / deconfinement phase transition 
for A/" = 4 supersymmetric Yang-Mills on the sphere. A slightly more complicated version 
was discussed in |13], where the transition between the low temperature and high temper- 
ature phase was shown explicitly by evaluating the actions of various compactifications of 
D4 branes. In the context of supergravity backgrounds dual to confining field theories the 
analog of the Hawking-Page phase transition was observed numerically in ^5j, similar results 
were reported in [16] . In every case the dominant thermal phase is determined by comparing 
the supergravity actions. 

The problem we have at hand, that is, the competition between the backgrounds where 
temporal boundary conditions are imposed along the r cycle and along the </> cycle seems to 
have the same flavor. However, a naive evaluation of the difference of actions, which in our 
case, as in the case of [13], boils down to an evaluation of the difference in volumes yields a 
very puzzling result. Namely, if we evaluate the difference in actions from the manifold with 
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temporal boundary conditions along r and along we obtain a result of the form: 



AV = Vr- V, 



oo 




f3o 



) 



(5.1) 



2 



L 



where /5o is the black hole inverse temperature. This result is puzzling because, since AS ~ 
AV, it seems to suggest that for high temperature the background (or the phase) that 
dominates is the one with temporal conditions along </>. Our intuition tells otherwise, we 
believe that the high temperature phase that dominates is the Schwarzschild black hole. We 
do not have a full resolution of this puzzle and hope to return to it in a separate publication. 
However, there are a few things that might go wrong with the above naive evaluation. For 
instance, the conical singularity along the cycle prevents us from evaluating the action 
simply as a volume since the curvature will develop singularities. 

We finish this section mentioning another idea that evolved largely from our previous 
work [32] and correspondence with T. Takayanagi. Namely, that the calculations presented 
in this paper, which fully involved the r-cycle, are to be understood as the Polyakov loop 
type of generalization of the standard entanglement entropy. This is a thought that we also 
plan to develop separately but that seems worth mentioning given the substantial difference 
between the phase structure observed in [32] with generically no phase transitions and the 
rich structure encountered here. 

6 Conclusions 

In this paper we have investigated the holographic proposal for the geometric entropy by 
computing the entanglement entropy for Schwarzschild black holes in global AdSp spaces. 
In the Euclidean version of the gravity geometry we have various choices of the subsystem 
A. Our calculations show, as already suggested in [31], that the geometric entropy we 
compute here is the analog of the Polyakov loop in field theory, whereas the calculations of 
the entanglement entropy in [32] are more like the analog of the calculations of the Wilson 
loop. 

For AdS^ we found a perfect agreement with the form of the geometric entropy computed 
for a two-dimensional CFT with periodic boundary conditions. We have also found (and 
believe that it is one remarkable result of our investigation) an universality of the phase 
diagrams for both AdS^ and AdS^. 

When interpreted from the field theory point of view our results fits many expectations 
about the structure of divergences and the various dependence of the entanglement entropy. 
In particular, we found an agreement with various expectations in 2 + 1 dimensions [37] and 
[36] . Similarly, our results fit nicely with expectation for 3 + 1 dimensional field theories as 
presented in [21 |T5] . However, due to the fact that we introduce various scales in the problem, 
our results have different N scaling than the expected from a conformal theory. We discussed 
this in the context of the dual field theories. We believe it would be interesting to pursue 
the implications of these scalings as potential phases that have not been considered before. 
In a sense this situation is similar to generalizations of the UV/IR relation in the context 
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of the AdS/CFT. It would be interesting to find, as in [SIIIH], precise descriptions on both 
sides of the correspondence of probes which yield different energy- in-the-boundary/distance- 
in-the-bulk relations. 

One interesting problem that we also plan to pursue is the following. From the field 
theory point of view our calculation corresponds to the geometric entropy for a field theory 
in a thermal state |\&). In much of the literature the state |\E') has been taken to be the 
ground state. What we found out in the present paper and in our previous explorations |32j . 
is that not only the state l^') is crucial but also that the region A plays a central role. It 
would be interesting to compute the geometric entropy at finite temperature directly in field 
theory. 
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A Piece-wise smooth minimal surface 

We want to determine the minimal piece-wise smooth surface that is connected to the horizon. 
Given the metric 

ds"^ = T{r)df + h{r)dr'^ + s(r, y')dx^ + atjir, y')dy'dy^ (A.l) 

with ro < r < oo, h{r Tq) ^ oo and r(r — rg) ^ 
We can parametrize the surface as: 

—1/2 < X < a s.t. Xa{ro) = a 

(A.2) 

b<x <l/2 s.t. Xb{ro) = b. 

Xa and Xb are determined by minimizing their area under the given boundary conditions. 
The minimal surface connecting Xa and Xb at ro is given by the segment on the horizon. We 
start by proving this. Parametrize the surface as 

t = constant, r = r{x) with b.c. r(a) = r{b) = Tq. (A. 3) 

The induced metric is 

ds'^ = [s(r, y') + h{ry^{x)]dx^ + atjdy'dy^ (A.4) 



= t-tr 



X — Xa{r) oo < r < ro, 
r — ro a < X < b 
X — Xb{r) ro < r < oo. 
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with area given as 

A = j^dx j dY^^s + (A.5) 

where dY ^/a is the volume element of the internal manifold. The problem of finding the 
minimal surface reduce to finding the equation of motion for a point particle r{x) with 
Lagrangian 

L = j dY^^s + /ir'2. (A.6) 
Since this Lagrangian does not explicitly depend on x, the quantity 

E = - i (A.7) 

is constant. Thus the e.o.m. is 

E = -l 4^^. (A.8) 

Since h{r — > tq) — > oo, then a nonzero value of E requires h{r{x — > a))r''^{x — > a) = 
h{r{x — s> 6))r'2(x ^ 6) = c where c is a constant, thus r'(a) = r'(6) = 0. So E is simply given 
by 

E = - [ siro,y'')V(^{ro,y')dY 
J y^s{ro,y') + c 

Since r(a) = r{h) = tq, there exist a Xq in the interval (a, 6) where the function r{x) hits 
a maximum value, r^. if is smooth, then r'(xo) = thus 

E = - I y/sirm,y')a{rm,y')dY. (A. 10) 



(A.ll) 



We then have 

s{ro,y')y/a{ro,y^)dY f ^ — 

We observe that 

s(ro, y')^a{r^,t)dY \ < J ^ s{r^,t)a{ro,y')dY if c > 



^s{ro,t) + c \> j ^s{rQ,t)a{rQ,t)dY if c < 
If the quantity 



(A.12) 



^/W)= / ^s{r,t)a{r,y^)dY (A.13) 



is increasing in r, then > tq exist only when c < 0. Since c cannot be negative, then r„ 
must be equal to tq. So we observe that the minimal surface is given by r(x) = tq. 
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Now we determine Xa and x^. In what follows, we drop the dependence of y in s{r,y^). 
As noted above, the y dependence is irrelevant. The induced metric under these constraints 
is 



e{r) + s\r) 



dXn 



dr 



dr 



dr'^ + aijdy^dy- 



dr^ + aijdy^dy 



(A.14) 



The area of the surface is 

A 



int 



ro 



h\r) + s\r) (^-^ 



1/2 



dr. 



(A.15) 



Vint is the volume of the internal manifold. We observe that the problem of finding the 
minimal area reduces to solving for the classical path of a point particle Xa- Since the 
Lagrangian does not explicitly depend on the position of the particle, then the Euler Lagrange 
equations implies 



d_ 

dr 



PL 
dx' 



where L = Vint\/h'^{r) + s^(r)(x'^)^ 



(A.16) 



where ' denotes derivative with respect to r. The equation of motion for Xa is given by 



/2( 2ta2 _ 2 



(A.17) 



where c is an integration constant to be determined by the boundary condition at r = tq. 
However we note that the r.h.s. tends to infinity as r ^ tq. This implies that the derivative 
of Xa blows up at the horizon whenever c ^ 0. The solution that minimize the action 
corresponds to then c = 0, i.e. 

Xa{r) = constant, (A. 18) 

and similarly for X},{r). To match the boundary conditions at the conformal boundary, we 
find that the only allowed discontinuous surface for the space-times with zero volume end 
space is: 

{X + 1/2 oo < r < ro, 
r-roa<x<b . (A.19) 

X — 1/2 ro<r<oo 

B Geometric entropy from AdS^i r{9) embedding 

We start by proving that AA changes sign. This quantity is given by, 

y^/m 



AA = 2pL 



Vy^fiy) -yUiy* 



1 My - (y* - yo) 



(B.i; 



25 



First consider the case when is large; then /(y) can be substitute for its asymptotic value 
of y"^. We thus obtain 



= 2pLy. 



l\dx-l 



(B.2) 



1 

= 2(3Ly^ [.4009... - 1] < 0. (B.3) 

So in the large limit, the difference in area is negative. 
Now consider the y^ — > yQ limit. It is clear that 

AA(y, = 2/o) = (B.4) 

since /(yo) = 0. Now we want to show that for y* = yo + e, AA > 0. For this we can deduce 
that the difference in area changes sign since it is positive for y^, near y^ and negative for 
large y^. We show this by proving that the derivative of A A with respect to y^, is positive at 
y^ and near yo. To facilitate taking derivatives, we rewrite the above expressions in terms of 
the coordinate x = —: 

V 

f{y) = -^{yl + yy-x^y^ + yo)) (B.5) 
y'fiy) = -M + y*x'-xHyo + yo)) (b.6) 

y'fiy) - ylf{y*) = - x){x\yl + - yl - yo) + x\yl + y.) + ylil + x)).(B.7) 

In this coordinate, AA becomes. 



AA = 2f3Ly^!}( ^ui^y,^-hi^yo)^ _ A 

yx'^V'i^\J{yl+y*-yt-yo)x''i+{yl+y*W+(.^+^)y'i J 
- 2(3L{y,-yo). (B.8) 

Taking derivative of AA with respect to y^, at fixed yo 



dAA ^ x2i/«(-2(x+l)yg+3(x2+x+l)(i;g+yo)y^+x2(yg+i;o))dx 



|3/2 



+ ^-2/3L^. (B.9) 

y* y* 

First, we observe that when y* = yo, 

= 3L fl x^yf,il+3y^)ix^+yl{l+x+x^))dx ^gj^ 

dy^ ^ V^\/{'i^-x){yox'^+yl{l+x+x'^))[yox'^+yl{l+x+x'2)\^/'^ ^ 
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We observe that the first term is positive and logarithmically divergent due to the factor 



in the integrand. This divergence comes from the fact that the term a/— {Vq + yo) + y*x'^ + vl 
in the denominator of the integrand in equation IB. 81 vanishes when y^, = y^ and x = \. Fur- 
thermore, one should note that all the other terms in the integrand are regular within the 
integration range except for the a/1 ~ x term; in fact the integrand is zero at x = 0. Thus 
when is close to y^, most of the integration will be coming from x near 1 region. To 
illustrate this picture better, we rewrite the integral as 

h{y^, yo,x)da^ 



VI - a: 



— X 

yo 



where h{y^:,yo,x) is regular for x G [0, 1]. Notice that we have factored the term that leads 
to the divergence at y^, = yQ since 

- {yl + yo) + y.x' + yl = - x)(yo'(4 + + x') + y^x'). (B.12) 

yo yo yo 

Written this way, it is clear that when is near yo, most of the contribution of the integral 
comes from a; ~ 1. Thus for y^, close to yo, we can write 

-PLh{l, 2/0, 2/o) In ( ^ - 1 ) + — - 2/?L > 0. (B.13) 

^y* 

since — is finite. 



dy* \yo J y* 



C A limiting value of l{y^): r^r) embedding in AdS^ 

In this note we demonstrate that 

Kv* ^yo) = l (c.i) 

for r(r) in AdS4. 

The main quantity is given as 

i( ) = - r y*v^ ^ (c 2) 

pJy. fiy)Vy'fiy)-yUiy*) ^' 

we first transform to a more convenient coordinate x = — where the various functions 
becomes, 

f{y) = -^{y'.+y*x'-x'{yl + yo)) (C.3) 

y'fiy) = -M + y*^'-^'(yo + yo)) (c.4) 



x^ 

y* 

x^ 
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y'fiy) - ylf{y.) = - x)ix'iy! + y.- yl - yo) + x\y', + y,) + y^l + x)).(C.5) 



From these equations, we obtain 



^ -^^ Vi^^/(x'Hy'i+y*-yo-yo)+x'^{yi+y*)+y'i{^+x)){yi+y*x2-x^{y^+yo))' 

This function has several characteristics that comphcate the analysis near ?/o- These are 

1. The pre-factor of the integral vanishes when y^, = yo. 

2. The denominator of the integrand have a zero when a; = 1. Since this zero comes under 
a square, it is not a serious problem. 

3. When = y^, the integrand is unbounded in the x = 1 limit. From numerical analysis, 
we know that as y^, approaches yo, this unboundedness is such that the integral diverges 
the right way to cancel the zero in the pre-factor. 

The key to obtaining an analytic understanding of / in the y^, yo limit, is to have a 
way of making case 3 precise. We can do this by factoring out the terms in the integral that 
are responsible for these pathologies. We have 



yl + y*-yo-yo = {y* - yo) {yl + y*yo + yl + i) (C.7) 

yl + y*x^ - x^{.yl + yo) = {.y*-yQx){yl + xy^yQ + x^yl + x^) (C.8) 



Now we can write /(y*) as 



l{y.) = l Hi^lMMy^y^dx (C.9) 
(3 Jo V 1 - X y*- yox 



where 



vix,y*,yo) 



X 



yWyo + y*yo + y* + 1 



{y^x^ + x2 + yoy^x + yl) a/ {-yl -yo + yl + y*) + {y^ + y^) + (x + l)yf 



Now we have isolated all the pathologies and introduced the function rj. This function is 
everywhere regular. This allows us to take limits in an organized way. Let 

e = ^^— ^. (C.IO) 

yo 

Now we perform a set of coordinate transformations on the integral equation /(y*) 

^ ^ 2 r]{x, y,, yo) y/y, - ^o ^^ ^ r]{x,y^,yo) 1 (C 11) 

p Jo ^Jl-x y^ - yox p Jo ^Jl-x - yox 

2 [T /S(l-z,,.,,o) 1 ^^^^^ 



/3\l yo Jo V ^ e + 2; 

28 



We now consider a coordinate change z = u"^ and the integral becomes 



Using that 



1 = ^ h^^-'^'^i-y'^du (C.13) 

Py/m Jo e + m2 



1 

arctan ( ^ ) , (C14) 



we are ready to take the e — >■ hmit in the expression for /. Now we let go to yo and e go 
to 0. We then obtain 

Ky* yo) = ^— — 77- (C.15) 



Now we have 



r/(l,2/o,2/o) = (C.16) 



This allows us to conclude that 



Ay.^yo) = \ (C.18) 

This sort of analytical result is very useful because it is used as a calibration for the 
numerics that we performed. In can clearly be seen from figure [5] that at the minimal value 
of which is y^^ the normalization of / is precisely as analytically obtained above. 
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